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This  paper  is  an  attempt  to  investigate  the  nozzle  expansion  characteristics  of  a  solar 
thruster  that  uses  ammonia  as  propellant.  For  this  purpose  we  have  developed  a  state  to 
state  kinetic  model  to  study  the  dissociation  of  ammonia  in  supersonic  nozzle  expansion.  The 
properties  of  the  thruster  depend  on  the  ammonia  dissociation  degree  in  the  reservoir,  which 
can  be  very  far  from  equilibrium  due  to  the  very  slow  dissociation  kinetics  at  the  thruster 
working  conditions  (<2000  K).  We  have  extended  the  calculation  for  higher  temperatures 
(>2500  K)  where  ammonia  dissociation  is  higher.  Due  to  the  presence  of  some  hydrogen 
atoms,  non-Boltzmann  distributions  have  been  observed. 
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Nomenclature 

nozzle  section 
Plank  constant 

constant  pressure  heat  capacity 
constant  volume  heat  capacity 
speed  of  light 
reaction  activation  energy 
electron  charge 

rate  coefficient  of  the  direct  reaction 

equilibrium  constant 

rate  coefficient  of  the  reverse  reaction 

reservoir  pressure 

universal  gas  constant 

reservoir  temperature 

vibrational  temperature 

vibrational  quantum  number 

Cv/R 

ammonia  dissociation  degree 
vibrational  level  energy 

Cp/Cv 
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me  =  first  order  spectroscopic  constant  to  calculate  vibrational  level  energies 

o)eXe  =  second  order  spectroscopic  constant  to  calculate  vibrational  level  energies 

meye  =  third  order  spectroscopic  constant  to  calculate  vibrational  level  energies 

coeZe  =  fourth  order  spectroscopic  constant  to  calculate  vibrational  level  energies 

I.  Introduction 

ecently,  solar  thermal  propulsion  (STP)  engines  have  been  investigated  as  a  promising  technology  to  enhance 
the  capabilities  of  micosatellites  maneuverability  and  enable  new  mission.  Solar  thermal  propulsion  takes 
advantage  of  concentrated  sunlight — typically  at  concentration  ratios  of  up  to  10000:1 — in  order  to  raise  a  heat 
exchanger  to  temperatures  of  several  thousands  degrees  K.  At  the  Surrey  Space  Centre1'2  (SCC),  a  small  scale  STP 
engine  has  been  designed  to  use  ammonia  as  propellant.  The  main  advantage  of  ammonia  as  propellant  is  that  it  is  a 
storable  liquid  with  low  molar  mass  and  easy  handling  characteristics  (and  does  not  exhibit  the  serious  toxicity  and 
flammability  risks  posed  by  hydrazine). 

To  achieve  improved  performance,  the  ammonia  gas  is  rapidly  heated  to  2000-2500  K,  which  should  result  in 
dissociation  to  N2  and  H2  molecules,  doubling  the  particle  density  and  decreasing  the  average  molecular  weight  of 
the  mixture.  This  improves  rocket  performance  by  raising  specific  impulse  (Isp),  which  varies  as  the  inverse  square 
root  of  molecular  weight.  The  problem  is  that  the  dissociation  kinetics  at  temperatures  of  interest  in  the  engine  are 
very  slow  and  therefore  the  complete  dissociation  of  ammonia  molecules  may  not  occur. 

In  this  paper  we  have  developed  a  kinetic  model  for  ammonia  mixtures  applied  to  the  nozzle  investigated  at 
SSC1’2.  It  is  worth  noting  that  we  have  introduced  state-to-state  kinetics  to  verify  the  presence  of  non-equilibrium 
vibrational  distributions  during  the  nozzle  expansion.  We  have  also  investigated  the  behaviors  of  the  engine  for 
higher  reservoir  temperatures  (3000K-5000K)  where  the  dissociation  of  ammonia  is  complete. 

II.  Model  Description 

A.  Nozzle  Geometry 


The  nozzle  used  in  the  thruster  has  a  conic  geometry  of  the  following  dimensions  (in  mm) 

throat  radius 

0.35 

exit  radius 

2.0 

exit-throat  distance 

6.2 

inlet  radius 

1.0 

inlet-throat  distance 

1.2 

These  data  have  been  provided  by  SSC 
researchers.  The  converging  nozzle  is  quite  short 
so  we  choose  to  include  part  of  the  reservoir  in  the 
calculation  through  a  connecting  cone  of  the 
following  dimensions  (in  mm) 
reservoir  radius  1 1 .25 
connecting  cone  length  2.0 
The  nozzle  is  presented  on  fig  1 . 

The  fluid  dynamic  model  of  the  nozzle  flow  is 
based  on  the  almost  one-dimensional  Euler 
equations  for  a  variable  section  duct.  The  Euler 
equations  are  solved,  for  the  stationary  state, 
together  with  the  perfect  gas  equation  of  state  and 
the  mass  continuity  for  each  species  including  the 
reaction  contributions.  The  numerical  scheme 
adopted  is  described  in2. 

B.  Chemical  Model 

To  build  the  chemical  kinetic  models,  we  need 
to  select  both  the  species  that  must  be  taken  into 
account  and  their  mutual  reactions.  In  the  initial 
demonstration  laboratory  experiments,  the 
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pressure  ranges  from  1  to  4  atm  and  the  temperature  is  less  than  1500  K.  For  operational  versions,  the  temperature 
could  reach  the  value  of  2500  K.  The  pressure  values  are  measured  directly,  but  the  gas  temperature  is  estimated 
from  some  measurements  performed  outside  the  reservoir.  This  results  in  a  large  uncertainty  on  the  actual  conditions 
operating  inside  the  reservoir.  For  these  temperatures  the  ammonia  decomposition  kinetics  are  quite  slow.  Surface 
dissociation,  neglected  in  this  approach,  should,  probably,  be  inserted  in  the  model. 

1.  Species 

We  assume  that  the  ammonia  mixture  consists  of  atoms  and  molecules  formed  in  the  NH3  decomposition.  We  ignore 
all  the  ions  because  the  temperatures  involved  are  quite  low.  We  introduce  the  vibrational  excitation  only  for  N2  and 
H2  molecules,  which  have  stable  vibrational  excited  levels.  For  all  the  other  molecules,  we  diregard  the  internal 
levels  as  their  number  is  so  high  that  the  computational  load  becomes  unaffordable.  A  multitemperature  model  could 
be  a  valid  approach  assuring  an  acceptable  computational  cost  nevertheless,  in  this  work,  we  consider  these  internal 
states  in  equilibrium  with  the  translational  degree  of  freedom  or  as  following  a  frozen  kinetics.  According  to  the 
internal  state  model  selected  we  assign  the  Cp  and  y  (y  =  Cp/Cv)  coefficients  to  each  molecule.  For  each  vibrational 
mode  or  rotational  axis,  in  equilibrium  with  the  translational  degree  of  freedom,  we  add  the  quantity  aR  into  the 
expressions  of  Cp  and  Cv,  where  R  represents  the  universal  gas  constant  and  the  a  parameter  is  set  to  1  for  a 
vibrational  mode  and  to  !4  for  a  rotational  axis.  This  approximation  is  valid  for  the  rigid  rotor  and  for  the  infinite 
harmonic  oscillator,  since  in  general  the  Cp  coefficient  depends  on  the  temperature.  For  atomic  species  no  internal 
degree  of  freedom  is  considered.  The  contribution  of  translational  degrees  to  the  molar  specific  heat  is  1.5R.  Table  1 
reports  the  species  accounted  and  their  relative  properties. 

Tab.  1:  species  inserted  and  their  physical  properties.  Rot  is  the  number  of  independent  rotational  axis,  Vib  the 
number  of  vibrational  modes  and  Ef  the  enegy  of  formation.  _ _ 


Species 

Mass  (AMU) 

Rot 

Vib 

Ef(eV) 

nh3 

17 

3 

6 

-0.403242 

n2h3 

31 

3 

9 

4.30994 

nh2 

16 

3 

3 

2.00285 

NH 

15 

2 

1 

3.90225 

n2 

28 

2 

1 

0.00000 

h2 

2 

2 

1 

0.00000 

N 

14 

0 

0 

4.87950 

H 

1 

0 

0 

2.23910 

For  H2  and  N2  we  estimate  the  energy  of  each  vibrational  level  (ev)  through  the  semi-empirical  formula,  function  of 
the  vibrational  quantum  number  and  of  some  spectroscopic  constants  expressed  in  cm"1: 


S 


V 


G)gxe 


+  0)ey  g 


+  (OgZg 


(1) 


The  H2  spectroscopic  constants  have  been  obtained  through  the  fits  of  the  vibrational  level  energies  extracted  from 
the  BKMP  potential  energy  surface4  while  for  N2  we  used  data  from  ref.  5.  Table  2  summarizes  those  values  for  the 
two  species. 


Table  2:  Spectroscopic  constants  (in  cm"1)  used  in  the  calculation  of  the  vibrational  level  energies. 


Parameter 

n2 

h2 

0)g 

2358.57 

4464.7 

«exe 

14.324 

146.93 

weYe 

-0.00226 

4.7108 

(OgZg 

0.0 

-0.22712 

Ediss 

9.7639 

4.4772 

V  max 

48 

15 
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2.  Processes 

We  implement  in  our  model  macroscopic  kinetics  for  NH,  depletion  based  on  the  following  processes: 


NH,+X  NH2+H+X 

(pl) 

NH,  +  X  ^  NH  +  h2  +  X 

(p2) 

NH3  +  H  NH2  +  h2 

(p3) 

NH  +  H  +  X  NH2  +  X 

(p4) 

H2  +  N  +  X  <H>  NH2  +  X 

(p5) 

NH2  +  H  NH  +  H2 

(p6) 

NH2  +  NH2  NH,  +  NH 

(p7) 

N  +  H  +  X  NH  +  X 

(p8) 

NH  +  H  H2  +  N 

(p9) 

NH  +  N  N2  +  H 

(plO) 

NH  +  NH  ^  H2  +  N2 

(pH) 

2H  +N2  H2  +N2 

(P 1 2) 

2H  +  H2  H2  +  H2 

(pl  3) 

2N  +N2  «-*■  N2  +N2 

(P 14) 

2N+N  ^  N2  +  N 

(p!5) 

nh3  +  nh2  n2h3  +  h2 

(P 1 6) 

N2H3  +  X  -*  NH2  +  NH  +  X 

(Pl7) 

N2H3  +  H  -*  2NH2 

(pl  8) 

N2H3  +  H  -»  NH3  +  NH 

(pl  9) 

where  X  is  a  generic  component.  The  rate  coefficients  for  all  these  processes  are  expressed  in  the  Arrhenius  form: 

K  =  K0T6e“Ea/T  (2) 

For  the  direct  reaction,  the  parameters  K0,  6  and  Ea  are  reported  in  Tab.  3.  The  rate  coefficients  of  these  processes 
have  been  taken  from  ref.  5-10.  In  this  model,  we  neglect  more  complex  molecules,  such  as  N2H4,  that  do  not  seem 
to  be  significant. 

Rates  for  the  reverse  processes  have  been  obtained  through  the  detailed  balance  principle 


—  - 

Keq  (3) 

where  Keq  is  the  equilibrium  constant.  Equilibrium  constants  have  been  calculated,  from  partition  functions, 
following  the  classical  statistical  thermodynamic  theory11.  Subsequently  they  have  been  fitted  by  the  equation 

T 

ln(K^)  =  K„  +  KpH  +Ktel 

\  1  /  (4) 

We  consider  only  the  equilibrium  constants  for  the  generation  of  a  new  species  from  standard  ones  like  N2  and  H2  in 
nitrogen-hydrogen  compounds  formation.  Therefore  all  the  formation  reactions  can  be  synthesized  as 


n 

2 


N,+ 


“H2-»N„Hm 


The  parameter  values  utilized  to  calculate  Keq,  for  the  species  listed  in  table  1  are  reported  in  tab.  4. 


(5) 
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Table  3:  Arrhenius  coefficients  for  the  reactions  inserted  in  the  model.  K0  has  the  dimension  of  cm3(n"1Vs,  where  n 
is  the  number  of  particles  involved  in  the  reaction,  Ea  is  expressed  in  K  and  6  is  a  dimensionless  parameter  (eq.  2). 


process 

K0 

6 

Ea 

Ref. 

pi 

1.53e-8 

0 

42400 

6 

p2 

1.00e-9 

0 

47000 

4 

p3 

3.20e-13 

0.67 

1720 

5 

p4 

1.00e-32 

0 

0 

7 

p5 

1.00e-34 

0 

0 

7 

p6 

2.30e-13 

0.67 

2160 

5 

P7 

6.64e-12 

0 

2800 

6 

P8 

l.OOe-33 

0 

0 

7 

P9 

1.70e-12 

0.68 

957 

5 

plO 

1.77e-l  1 

0 

0 

5 

pll 

6.60e-13 

0.55 

957 

5 

pl2 

1.34e-31 

-0.60 

0 

5 

pl3 

2.68e-31 

-0.60 

0 

5 

pl4 

7.44e-32 

-0.50 

0 

5 

pl5 

3.31e-27 

-1.50 

0 

5 

pl6 

1.33e-12 

0.5 

10850 

4 

Pl7 

1.70e-8 

0. 

21000 

4 

p  1 8 

2.60e-12 

0 

0 

4 

pl9 

1.70e-13 

0 

0 

4 

Table  4:  Coefficients  used  to  calculate  the  formation  equilibrium  constants  when  the  pressures  are  expressed  in  Pa. 


Species 

KP 

qP 

Ke 

9e 

nh3 

51.8293 

-13.1631 

0.960208 

-7.97367 

628.095 

nh2 

19.6666 

63.4542 

1.01863 

13.5452 

311.287 

NH 

-5.14978 

85.2911 

0.998065 

-2.21874 

423.673 

N 

28.3772 

-114.298 

0.997620 

-2.14877 

4371.66 

H 

26.1506 

-52.7759 

0.996362 

-2.71911 

2157.34 

From  these  formation  equilibrium  constants  we  can  calculate  the  equilibrium  constants  relative  to  all  the  reactions 
involving  the  species  in  the  table. 

The  kinetic  model  has  been  improved  including  state-to-state  kinetics  for  the  diatomic  species  N2  and  H2, 
substituting  processes  P12-P15  with  the  state-selective  dissociation  and  adding  vibrational  relaxation  processes12"14 


N2(v)  +  N2(w)  «  N2(v-1)  +  N2(w+1) 

(p20) 

N2(v)  +  N2(w)  «-*  N2(v-1)  +  N2(w) 

(p21) 

N2(v)  +  N  **  N2(v-Av)  +  N 

(p22) 

N2(v)  +  N2(w)  «-»  N2(v-1)  +  2N 

(p23) 

N2(v)  +  N2(w)  2N  +  N2(w) 

(p24) 

N2(v)+N^  2N+N 

(p25) 

H2(v)  +  H2(w)  **  H2(v-1)  +  H2(w+1) 

(p26) 

H2(v)  +  H2(w)  **  H2(v-1)  +  H2(w) 

(p27) 

H2(v)  +  H  «•  H2(v-Av)  +  H 

(p28) 

H2(v)  +  H2(w)  ■*-»  H2(v-1)  +  2H 

(p29) 

H2(v)  +  H2(w)  **  2H  +  H2(w) 

(p30) 

H2(v)  +  H  «•  2H  +  H 

(p31) 

All  the  other  processes  involve  only  the  ground  state  of  such  molecules, 
relaxation  of  nitrogen  molecues  due  to  hydrogen5'8: 

We  have  also  included  vibrational 

N2(v)  +  H2(w)  «-»  N2(v-1)  +  H2(w) 

(p32) 

N2(v)  +  H  N2(v-Av)  +  H 

(p33) 

N2(v)  +  H2(w)  **  N2(v-1)  +  H2(w+1) 

(p34) 
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III.  Results 
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inlet 
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C.  Low  reservoir  temperature  (T<2500  K) 

In  this  paragraph  we  want  to  investigate  the  kinetic  behavior  of  the  ammonia  mixture  in  cases  close  to  the 
working  conditions  of  the  solar  thruster  studied  at  the  Surrey  Space  Centre1'2.  The  pressure  range  is  between  1  and  4 
Bar  and  the  temperature  is  <  1500  K.  The  heater  temperature  is  not  well  known  because  it  has  been  mesured1'2 
outside  the  chamber.  We  have  extended  our  calculations  to  temperature  values  below  2500  K,  which  may  be 

obtained  in  orbital  conditions. 

We  investigate  four  different 
pressure  values  in  the  range  1-4  Bar. 
In  figure  2  we  report  the 
concentration  profiles  along  the 
nozzle  obtained  for  reservoir 
conditions  of  P=4  Bar  and  T=1500  K, 
and  introducing  the  vibrational 
kinetics  of  N2  and  H2. 

It  is  interesting  to  observe  that  the 
flow  is  practically  frozen  along  the 
whole  nozzle  except  for  some  minor 
species  such  as  H  and  NH2  that 
completely  disappear  at  the  nozzle 
exit.  The  atomic  nitrogen  molar 
fraction  is  lower  than  10"10. 

Figures  3  and  4  describe  the 
vibrational  distribution  of  hydrogen 
and  nitrogen  for  the  conditions  of 
figure  2.  We  can  observe  that  the 
nitrogen  distributions  are  frozen.  On 
the  contrary,  while  H2  distributions, 
for  v<3,  are  almost  frozen,  their  tails 
are  strongly  depleted,  resulting  in 
non-Boltzmann  distributions.  The 
different  behavior  of  the  two  species 
is  due  to  the  different  percentage  of 
nitrogen  and  hydrogen  atoms; 
hydrogen  atom  concentration  is  high 
enough  so  that  the  VT  collisions 

involving  atoms  are  frequent  enough  to  cool  the  distribution  tails.  We  observe  that  the  Mach  number  is  practically 
independent  of  the  initial  temperature  and  pressure,  obtaining  at  the  nozzle  exit  a  value  close  to  Mach  5.4.  We  report 
the  pressure  along  the  nozzle  in  figure  5.  It  is  quite  interesting  to  note  that  the  exit  pressure  is  proportional  to  the 
inlet  pressure.  These  two  results  occur  as  the  chemistry  weakly  affects  the  macroscopic  quantities  during  the 
expansion. 

Table  5  reports  the  mass  flow  rate  and  the  thrust  for  different  inlet  pressures  at  T0=1500  K  or  T0=2500  K.  The 
agreement  with  the  experimental  results  reported  in  table  6  is  satisfactory,  as  the  nozzle  reservoir  conditions  (mainly 
the  temperature)  are  not  well  known  experimentally. 

In  this  temperature  range  the  kinetics  of  ammonia  presents  large  uncertainties5"8,  even  though  it  has  been  widely 
studied.  The  main  problem  is  that  the  slow  kinetics  prevents  us  from  estimating  the  kinetic  constant  precisely. 
Experimental  works1'2  on  ammonia  thrusters  assume  a  complete  decomposition  of  ammonia  in  the  reservoir,  but  the 
kinetic  scheme  proposed  does  not  permit  to  reach  the  equilibrium  condition  in  an  acceptable  time  for  T<2000  K. 
Under  these  conditions  the  most  important  problem  is  that  it  is  possible  that  in  the  inlet  conditions  there  is  not 
equilibrium  as  usually  happens  in  nozzle  reservoirs.  Because  of  this  we  also  try  as  initial  condition  a  pure  ammonia 
gas. 

Dissociation  inside  the  reservoir  and  the  nozzle  is  very  poor,  so  that  again  in  this  case  the  flow  is  frozen.  As  a 
consequence,  the  flow  properties  depend  only  on  the  Cp/Cv  ratio  and  on  the  mean  molar  mass,  so  the  exit  properties 
are  completely  different  form  the  equilibrium  case.  As  an  example  in  figure  6  we  compare  the  Mach  numbers 
computed  starting  from  an  equilibrium  composition  (Eq.)  and  from  pure  ammonia  (NoEq.).  It  is  not  possible  to 
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Figure  2:  Molar  fractions  of  mixture  species  along  the  nozzle.  We 
assume  as  reservoir  conditions  Po=4  Bar  and  T0=1500  K  and  start  from 
an  equilibrium  inlet  composition. 
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really  know  what  composition  prevails  inside  the  reservoir  for  T<2000K,  moreover  some  impurites  or  wall  catalytic 
effects  can  accelerate  the  kinetics,  bringing  the  system  to  equilibrium  faster  than  predicted  by  our  kinetic  model. 


Fig.  3:  Hydrogen  vibrational  distributions  at  different 
nozzle  positions  in  the  same  conditions  as  in  fig.  2. 


Figure  4:  Nitrogen  vibrational  distributions  at  different 
nozzle  positions  in  the  same  conditions  as  in  fig.  2. 


T  =1500 K 

o 


inlet  throat 


x(m) 

Fig.  5:  Pressure  profile  along  the  nozzle  as  a  function 
of  the  inlet  pressure  for  T0=1500  K  and  starting  from 
an  inlet  equilibrium  composition. 


Table  5:  Computed  mass  flow  rate  and  thrust  of 
the  nozzle  as  a  function  of  the  inlet  pressure.  The 
inlet  temperature  is  T0=1500  K  except  when 
marked  with  the  symbol  assigned  to  T0=2500 
K.  At  the  inlet  we  suppose  an  equilibrium 
composition. 


Po(Bar) 

mass  flow  rate  (g/s) 

thrust  (mN) 

1 

2.1708e-2 

64.337 

2 

4.3404e-2 

128.68 

3 

6.5077e-2 

192.97 

4 

8.6784e-2 

257.38 

4* 

6.5954e-2 

256.74 

*  T=2500 K 


Table.  6:  Mass  flow  rate  and  thrust  of  the  nozzle 
determined  experimentally  in  different  test  cases. 


mass  flow  rate 
(g/s) 

thrust  (mN) 

.140 

219 

.147 

223 

.117 

183 

.121 

190 
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T  =1500 K 

o 


x(m) 


Figure  6:  Mach  number  profile  compared  for  inlet 
equilibrium  (Eq)  and  pure  ammonia  (noEq)  case 
(pressure  and  temperature  conditions  as  in  lig  2). 


Table  7:  Computed  mass  flow  rate  and  thrust  as  a 
function  of  the  inlet  dissociation  degree  (aD)  of  the 
ammonia.  Inlet  pressure  and  temperature  are  Po=4  Bar  and 
T0=1500  K. 


aD 

mass  flow  rate 
(g/s) 

thrust  (mN) 

Isp 

0.0 

1.1 189e-l 

284.13 

259 

0.2 

1.0437e-l 

280.35 

274 

0.4 

9.7946e-2 

275.67 

287 

0.8 

8.9483e-2 

264.37 

296 

1.0 

8.5450e-2 

253.33 

302 

We  try  the  same  comparison  at  T0=2500  K.  In  this  case, 
the  mixture  reaches  the  equilibrium  composition  in  the 
reservoir.  Nevertheless,  the  Mach  number  profdes  in  the 
equilibrium  and  non  equilibrium  case  are  again  different 
(see  fig.  7).  This  is  the  demonstration  that  the  history  of  the 
gas  in  the  reservoir  is  important. 


T  =2500  K 
0 


Figure  7:  Same  as  in  fig  6  with  reservoir 
temperature  T0=2500  K  and  pressure  Po=4  Bar. 


T  =1500  K 

0 


x(m) 


Figure  8:  Reduced  temperature  profile 

compared  for  inlet  equilibrium  (Eq)  or  pure 
ammonia  (noEq)  case  (pressure  and 
temperature  conditions  as  in  fig  2). 


Effects  of  initial  conditions  propagate  also  on  the  temperature  profiles.  In  figure  8  we  show  the  normalized 
temperature  profile  calculated  starting  from  an  equilibrium  condition  (Eq.)  and  from  an  undissociated  ammonia 
(NoEq.)  for  a  reservoir  temperature  of  1500  K  and  a  pressure  of  4  Bar.  We  observe  that  the  exit  temperature  is  much 
higher  in  the  case  of  the  equilibrium  conditions.  This  behaviour  could  be  due  to  the  lower  mean  molar  mass 
associated  with  the  mixture  in  equilibrium  composition  inside  the  reservoir. 

Different  behaviours  occur  when  the  reservoir  temperature  reaches  2500  K  (see  fig.  9).  In  fact,  at  this 
temperature,  the  reactions  in  the  inlet  are  faster  than  in  the  flow,  therefore  one  can  observe  an  increase  of  the 
temperature  due  to  the  mixture  dissociation.  Even  if  the  ammonia  energy  is  lower  than  that  of  nitrogen  and  hydrogen 
molecules,  the  compression  due  to  the  double  particle  number  density  produces  an  increase  of  the  temperature.  At 
the  nozzle  exit,  the  differences  are  lower  than  in  the  case  of  an  inlet  temperature  of  1500  K. 
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T  =2500  K. 

o 


We  find  similar  behaviours  for  the  relaxation  of  the  vibrational  temperature,  calculated  as  the  temperature  of  the 
first  excited  level  (figure  10).  Starting  from  non  equilibrium  conditions  we  can  observe  an  increase  of  the  vibrational 

temperatures  following  the  gas  temperature.  At  the  nozzle 
exit  the  vibrational  temperature  freezes  at  a  temperature 
higher  than  2300  K  just  after  the  throat.  N2  vibrational 
temperature  is  practically  frozen. 

As  a  consequence  of  the  previous  results,  we  can  state 
that  there  are  some  uncertainties  about  the  inlet  composition. 
To  demonstrate  the  sensitivity  we  have  reported  in  table  7  the 
calculated  mass  flow  rate  and  thrust  for  different  dissociation 
degrees  “aD”of  ammonia  for  the  process  NH3  =  1.5  H2  +  0.5 
N2  defined  as 


an  =  1- 


[nh3; 

1 

[nh 

3] 

+.5 

N 

Fig.  9:  Reduced  temperature  profile  compared  for 
inlet  equilibrium  (Eq)  or  pure  ammonia  (noEq) 
case,  for  a  reservoir  temperature  T0  =  2500K 
(pressure  and  temperature  conditions  as  in  tig  2). 


that  ranges  from  0  (pure  ammonia)  to  1  (N2  and  H2  mixture). 
We  note  that  both  quantities  decrease  when  the  dissociation 
degree  increases,  but  the  thrust  changes  more  slowly  than  the 
mass  flow  rate,  so  the  specific  impulse  will  be  higher  for  the 
dissociated  case. 


T  =2500  K  P  =4  Bar 


x(m) 

Figure  11:  Molar  fraction  profiles  calculated  with 
different  kinetic  models. 


x(ra) 

Figure  10:  Vibrational  temperature  (Tv)  profile  of  High  reservoir  temperature  (T>2500  K) 

N2  and  H2  in  the  same  conditions  as  in  figure  9.  *n  l^e  Previ°us  paragraph  we  have  studied  the 

behavior  of  the  ammonia  solar  thruster  in  the  working 
conditions  of  the  thruster  investigated  experimentally  in 
refs.  1,2.  In  this  paragraph  we  want  to  understand  the  behaviors  of  the  thruster  for  higher  temperature  ranges  (3000 
K  <  T0  <5  000  K).  For  temperature  T0>2500  the  ammonia  dissociation  is  very  fast. 

Moreover  in  the  considered  temperature  range  the  molar  fraction  of  atomic  hydrogen  is  not  negligible.  In  this 
conditions,  due  to  the  recombination  of  hydrogen  atoms  and  to  atom-molecule  VT  collisions,  non  equilibrium 
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vibrational  distributions  can  be  observed  and  non-Arrhenius  chemical  rates  are  obtained12,  inducing  anomalous 
behaviors  in  the  kinetics. 

Here  we  want  to  focus  on  the  differences  between  the  results  obtained  with  two  different  models: 


a  -  macroscopic  kinetics 
b  -  state  to  state  kinetics. 


x(m) 

Figure  12:  Velocity  profiles  calculated  with  different 
kinetic  models. 


In  fig.  11  we  have  compared  the  molar  fractions  of 
some  species  calculated  with  the  two  models.  It  is  possible 
to  observe  that  H2  molar  fraction  increases  when  the  state 
to  state  model  is  considered,  while  the  macroscopic  model 
predicts  a  frozen  concentration.  On  the  contrary,  NH2  and 
NH  molar  fractions  show  smaller  variations  when  the  state 
to  state  approach  is  used. 

This  behavior  affects  the  macroscopic  quantities  such 
as  the  gas  speed,  the  temperature  and  the  thrust.  In  fig.  12 
we  have  compared  the  axial  speed  profile  in  the  two 
kinetic  models  at  two  different  reservoir  temperatures. 
There  are  appreciable  differences  at  the  nozzle  exit 
especially  for  reservoir  temperature  T0=5000  K.  This 
effects  is  the  consequence  of  the  stronger  recombination  of 
hydrogen  atoms  in  the  H2  molecules  as  can  be  observed  in 
figure  11.  Similar  behaviors  can  be  observed  in  the  gas 
temperature  profile  (fig.  13)  and  in  the  thrust  (fig.  14).  An 
important  result  is  that  the  calculated  thrust  (and  therefo 
vibrational  kinetic  model  is  included  in  the  calculation. 


x(m) 

Figure  13:  Temperature  profiles  calculated  with 
different  kinetic  models. 


x(tn) 

Figure  14:  Thrust  profiles  calculated  with  different 
kinetic  models. 

:  the  performances)  is  higher  when  the  state  to  state 


IV.  Conclusion 

In  this  paper  we  have  investigated  the  chemical  kinetic  aspects  of  an  ammonia  mixture  for  application  to  solar 
thermal  thrusters.  We  have  obtained  results  that  contradict  the  common  assumption  of  considering  inlet  conditions 
in  equilibrium.  If  the  kinetics  of  the  mixture  is  very  slow,  as  in  the  case  of  ammonia,  the  gas  composition  can  be  far 
from  equilibrium.  At  low  reservoir  temperature  (T0<2500  K),  that  is  the  working  condition  for  the  thruster  in  the 
laboratory,  the  kinetics  of  internal  states  is  negligible,  and  the  macroscopic  and  state  to  state  models  gives  the  some 
results,  mainly  because  the  flow  is  almost  frozen.  On  the  other  end,  at  higher  reservoir  temperature  (T0>3000  K), 
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when  the  hydrogen  dissociation  is  important,  the  vibrational  kinetics  becomes  very  important.  In  these  conditions 
the  macroscopic  and  the  state  to  state  model  predict  different  results,  mainly  due  to  the  different  atomic 
recombination. 

We  must  emphasize  that  this  study  is  preliminary.  Some  properties  of  the  mixtures  have  been  neglected,  such  as 
the  vibrational  energy  of  the  NHX  molecules  and  the  formation  of  NH  molecules  by  vibrationally  excited  N2  and  H2. 
The  model  can  be  improved  including  a  multitemperature  approach  to  account  for  the  vibration  of  polyatomic 
molecules  and  by  adding  new  state  to  state  reactions  between  N2  and  H2. 
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